home *** CD-ROM | disk | FTP | other *** search
/ Info-Mac 4 / Info_Mac IV CD-ROM (Pacific HiTech Inc.)(August 1994).iso / Science / MAL ©P.f.Howden 1⁄1⁄89 / DIANEQNS.TXT < prev    next >
Text File  |  1989-02-23  |  3KB  |  76 lines
  1. 'TEST DIFFERENTIAL EQUATIONS to <MERGE> with DIANALOG by PASTING.
  2. 'Remove ' from each required line(esp.line 10).
  3.  
  4. '10 X(1,2)=2*X(2,0)-5*X(1,0)
  5. 11 X(2,2)=2*X(1,0)-2*X(2,0)
  6. 12 X(3,5)=X(3,1)+1
  7. 21 RETURN:'SPRINGS:2 2nd order lin EQNS+5th order.SEE TEXT for TIME PLOT.
  8.  
  9. '10 X(1,2)=SQR((2+X(2,0)*X(1,1))/8)
  10. 11 X(2,1)=SQR(1+X(1,1)^2)
  11. 20 RETURN:'CATENARY:2nd & 1st order EQNS.INPUTS=
  12. '2,2,0,.1,2,2,0,1,0,0,2,RET,0,2,0,-2,35,1
  13.  
  14. '10 X(1,2)=-X(1,0)
  15. 11 X(2,2)=-X(2,0)
  16. 21 RETURN:'CIRCLE:2 2nd order lin EQNS.INPUTS=
  17. ' 2 ,2,0,.1,2,0,1,2,1,0,0,2,CONT,40,10,0,10,0,1
  18.  
  19. '10 X(1,2)=1+X(2,1)
  20. 11 X(2,2)=-X(1,1)
  21. 21 RETURN:'ELECTRON:2 2nd order lin EQN.INPUTS=
  22. '2,2,0,.1,2,0,4,2,0,0,0,2,RET,40,5,0,1,35,1
  23.  
  24. '10 V1=SIN(100*3.14159*T):IF V1<0 THEN V1=0
  25. 11 X(1,2)=500000!*V1-102*X(1,1)-2200*X(1,0)
  26. 21 RETURN:'REPLACE DIANALOG LINES WITH CORRESPONDING:
  27. '170 FOR J=1 TO Q+1:PRINT"X";J;" SCALE MULTIPLIER M(";J;")=";:INPUT M(J)
  28. '590 X(2,0)=V1:FOR J=0 TO Q+1:V(J)=INT(V+M(J)*X(J,0)+Y(J)+.5):NEXT J:REM CAN PLOT OTHER VARIABLES
  29. '640 FOR J=0 TO Q+1:IF J=0 AND N0=10 THEN N0=0:GOTO 690
  30. '750 VY=275-K6*10-(YO-YM*T)*7
  31. 'ELECT FILTER:1 2nd order EQN:1,2,0,.0005,0,0,0,2,RET,0,.5,0,400,35,4
  32.  
  33. '10 X(1,2)=X(1,0)*X(2,1)^2-2/X(1,0)^2
  34. 11 X(2,2)=-2*X(1,1)*X(2,1)/X(1,0)
  35. 21 RETURN:'Replace DIANALOG lines with axes transform lines:
  36. '740 XA=X(1,0)*COS(X(2,0)):VX=(V+.5+XO+XM*XA)*7+1
  37. '750 YA=X(1,0)*SIN(X(2,0)):VY=275-K6*10-(YO+YM*YA)*7
  38. 'ORBIT:2 2nd order EQNS:2,2,0,.2,2,3(4.5->escape ),0,2,0,.3333333(.2->escape)
  39. '0,2,RET,10,4(2->escape),40,2,20,3  :Formula: X(1,0)=9/(2+COS X(2,0)).
  40.  
  41. '10 X(1,2)=K(0)*(1-X(1,0)*X(1,0))*X(1,1)-X(1,0)
  42. 21 RETURN:'Replace DIANALOG line 750 with line:
  43. 750 VY=275-K6*10-(YO+YM*X(1,1))*7
  44. '2nd Order VAN DER POL Eqn :1,2,0,.01,1,1,2,2,30,5,0,5,25,2
  45.  
  46. '10 X(1,2)=2*X(2,0)-5*X(1,0)
  47. 11 X(2,2)=2*X(1,0)-2*X(2,0)
  48. 21 RETURN:'SPRINGS:2 2nd order  lin EQNS.
  49. 'INPUTS=2,2,0,.1,2,5,0,2,0,0,0,2,RET,40,4,0,4,20,1
  50.  
  51. '10 K(1)=26: X(1,2)=-K(0)*X(1,1)-13*X(1,0)+K(1)
  52. 21 RETURN:'Replace DIANALOG line 750 with line:
  53. 750 VY=275-K6*10-(YO-YM*T)*7
  54. '2nd Order Damped Spring Mass Eqn :Inputs=1,2,.025,0,0,1,2,RET,0,10,0,8,35,1
  55.  
  56. 'SINGLE BURST DRIVING FUNCTIONS FOR DIANALOG [see text]:Change variables (other than T) so as not to clash in DIANALOG.
  57. '8 T0=4:T1=6
  58. '9 INPUT "T=";T
  59. '10 IF T<T0 OR T>T1 THEN Y=0:GOTO 12
  60. '11 Y=1
  61. '12 F=Y*(T+1-T0):PRINT "F=";F:'DELETE THIS PRINT "F" WHEN SOLVING EQNS.
  62. '14  PUT D.Es HERE ON LINES 14-19
  63. '19 GOTO 9:'  DELETE THIS LINE WHEN SOLVING EQNS.
  64.  
  65. 'REPETITIVE BURST DRIVING FUNCTIONS FOR DIANALOG etc.PRESET Y=0,1 IF INITIAL FUNCTION TO BE ON,OFF [see text]:Change variables (other than T) so as not to clash in DIANALOG.
  66. '8 T0=3:T1=6:T2=1:N=1
  67. '9 INPUT "T=";T
  68. '10 IF T>=N*T1-T2 THEN N=N+1:Y=1
  69. '11 IF T>T0-T2+(N-1)*T1 THEN Y=0:GOTO 13
  70. '12 Y=1
  71. '13 F=Y*(T+1+T2-(N-1)*T1):PRINT "F=";F:'DELETE THIS PRINT "F" WHEN SOLVING.
  72. '14 PUT D.Es HERE ON LINES 14-19
  73. '19 GOTO 9:, DELETE THIS LINE WHEN SOLVING EQNS.
  74.                 MACINTOSH LISTING
  75.  
  76.